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We demonstrate that the scatter in the luminosity relations of astrophysical objects can be used 
to search for axion-like-particles (ALPs). This analysis is applied to observations of active galactic 
nuclei, where we find evidence highly suggestive of the existence of a very light ALP. 
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Precision observations of the recent Universe are an 
invaluable test-bed for 'new physics' such as the exis- 
tence of new, low mass and/or weakly interacting parti- 
cles. Such a particle, the axion, was proposed in 1977 
by Pecci and Quinn to solve the strong CP problem pQ. 
Since then, other models which also feature (very) light, 
neutral spin zero, axion-like-particles (ALPs) have been 
proposed e.g. HH2H5- ^ ecen ^yi analyses of starlight po- 
larization [5], and the distribution and spectrum of high 
energy cosmic rays [5] have provided tentative evidence 
for ALPs. 

In this Letter, we present a new method for studying 
ALPs using astrophysical X/7-ray luminosity relations, 
which, when applied to observations of Active Galactic 
Nuclei (AGN) , provides the strongest evidence yet for the 
existence of ALPs. 

ALPs are scalar or pseudo-scalar particles, <f>, which 
couple to photons, A^, via the terms (cj)/M)F flu F tJ/l ' 
or ((j)/M)e fl u pcr F fJ -"F f " 7 in the Lagrangian respectively; 
Ffj, v = 2dh l A v \. We define m^, to be the ALP mass, and 
g 77 = 1/M is the coupling between ALPs and photons. 

In the presence of a background magnetic field, B, 
ALPs mix with photons. The probability that an ALP 
converts into a photon whilst travelling through a coher- 
ent magnetic domain of length L is [7] 
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eB 2 /M 2 ; h, c = 1. 



u 2 Y = Ana cm n c /m c 



is 



the plasma frequency; n c is the electron number density, 
m e the electron mass and a cm w 1/137 is the electro- 
magnetic fine structure constant, e = +1 for scalars and 



for pseudo-scalars; generally B 2 /M 2 \\m ( j 1 — w 2 j|| <C 1. 
The total flux (in ALPs and photons) is conserved by the 
mixing process, however photon number is not. 

ALP- photon mixing is constrained by a number of lab- 
oratory experiments (see Ref. [4] and references therein) 
but the tightest constraints come from the astrophysical 
consequences of ALPs (see Ref. [5] ) . 



We are concerned with very light ALPs, <C 
10~ 12 eV. For such masses, observations of the super- 
nova SN1987A limit g n = 10 u GeV/M < 1 @] for 
pseudo-scalars, whilst limits on new long ranges forces 
require gn < 1CP 16 for scalars. However, if ALPs are 
chameleonic, these constraints do not apply j3j |5], and 
the best constraint comes from the structure of starlight 
polarization: gn < 100 [5]. 

As photon number is not conserved ALP-photon mix- 
ing in the magnetic fields of galaxies and galaxy clusters 
alters the observed luminosity of astrophysical objects. 
In this Letter, we consider light that has travelled through 
N 3> 1 randomly oriented magnetic regions. This is 
true of light from many astrophysical sources, particu- 
larly those in galaxy clusters. We focus on the limit in 
which the mixing is strong, NP^^^ 3> 1, and frequency 
independent, iVA < tt/2. 

In this strong mixing limit, little or no circular polar- 
ization is produced [S] and power is spread with equal 
probability between the two photon polarization states 
71 and 72 and <f>. We take /( tot ) 1 / 2 u = (71,72, 4>) T to 
be real, then since ||u|| 2 = 1 is conserved u describes a 
point on a sphere, and any point is as likely as any other. 
Thus, after strong mixing, the normalized final state u is 
a random vector with: 
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where K, ~ U[— 1, 1), and the final flux in the photon 



field is 7 f (7) = (1 - K 2 )L ( - tot \ The initial photon flux is 

Ig 7 ^ and we assume L^ s» /( tot ). A state labelled by a 
real vector u is sufficient to describe any mixture of <f> 
with fully linearly polarized light. We extend this result 
to light with partial or no initial linear polarization by 
noting that any such state, with L^ m L^ tot \ can be 
written as a superposition of two real u state vectors [5]. 
Where pq is the initial degree of linear polarization, the 
final photon flux after strong mixing is: 
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where K 4 ~ U(-l,l). We define C = I^/I^l Av- 
eraging C over many different light paths through a 
large number of randomly oriented magnetic regions gives 
C = 2/3. This average reduction in the apparent lumi- 
nosity of astrophysical sources is well known [3] , however 
attempts to use it to constrain ALPs have, so far, been 
unsuccessful since the intrinsic luminosity of the astro- 
physical objects is not known with sufficient accuracy. It 
is rarely appreciated that C = 2/3 only when averaged 
over many different light paths (e.g. light from many 
different objects). For a given light path through the 
magnetic regions, there is significant scatter and skew 
of C about its mean of 2/3 [5]. The probability that 
C € [c — dc, c + dc] is /c(c) dc where: 
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where a = (1 — po) /(l + po) and c± — min(c, (1 ± po) /2. 

The central idea of this Letter is that the scatter in em- 
pirically established luminosity relations can constrain, 
detect or rule out strong mixing between ALPs and pho- 
tons. Provided lo is high enough, strong and (almost) 
frequency independent mixing can be caused by the mag- 
netic fields of galaxy clusters [5] , the existence of which 
is well established [TU]. When such mixing occurs, there 
will be a contribution to the scatter of the observed lu- 
minosities as described by Eq. (|3]). Observations im- 
ply that most clusters contain magnetic fields of strength 
B«l - which are generally coherent over L ~ kpc 

scales, although L may be as large as 10— 100 kpc in some 
cases [10] . Typical electron densities in the diffuse plasma 
in the intra-cluster medium (ICM) are n e ~ 10~ 3 cm~ 3 



and hence 
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: eV, and m 2 « w£ if 9ll ~ 0(1) 



This last condition requires the ALP coupling to be very 
close to the upper bound from supernovae or for the ALP 
to be chameleonic in which case gn > 1 is allowed. 

If light travels a typical distance of 0.1 — IMpc through 
the ICM then N ~ 100 — 1000 magnetic domains will have 
been traversed. Strong mixing requires iVP 7w ^ 3> 1, 
so VN(BL/2M) > 1 and u > Mu^^NB. With 
B 5 = B/fyG, L 2 = L/2kpc and g n = 10 n GeV/M we 
have BL/2M « 0.l5g n B 5 L 2 . Typically B 5 ,L 2 ~ 0(1) 
and for N ~ 100 - 1000 and we need u > 16 - 51 eV. 
Frequency independent mixing, NA < 7r/2, i.e. requires 
lo > N(2tt/lo 2 1 L) = 3 — 30keV. Numerical simulations 
show that the frequency independent limit is still approx- 
imately valid for frequecies that are a factor of 10 smaller 
i.e. u > 0.3 - 3keV. When NP 1 ^ 4> > 1, iVA > 1 the 
measured luminosity is attenuated by a factor of 2 /3 with 
relatively little scatter. 

For the remainder of this Letter, we assume that 
lo > 2keV light, i.e. X/7 rays, which originated in or 



passed through a galaxy cluster, has undergone strong 
and (almost) frequency independent ALP mixing, requir- 
ing gn > 0.1-0.3 andm^, < 10~ 12 eV. For lo < 0.5 keV, 
mixing is highly frequency dependent and can be either 
weak (so C ~ 1) or strong. In the latter case the lumi- 
nosity is reduced by a factor ss 2/3. 

We require that the X/7 sources are compact i.e. their 
size, R, is < L ~ few kpc. Diffuse light, such as the 
X-ray light from galaxy clusters (R ~ 10 2 — 10 3 kpc) is 
not suitable for our analysis as it will have traced many 
different paths through the magnetic fields of the cluster, 
and the effects of mixing will be averaged over all of these 
paths. For such sources, the only effect of strong mixing 
is the 2/3 suppression of the total luminosity. 

For a number of classes of compact astrophysical ob- 
jects, correlations between the X or 7 ray luminosity or 
radiated energy and some feature of their light-curve (e.g. 
peak energy) or the object's luminosity at a lower fre- 
quency have been empirical established. We let Yi label 
the X or 7 ray luminosity / total energy and Xi label 
the light-curve feature or lower energy luminosity with 
which it is correlated. The relations between Yi and Xi 
take the form: 



log 10 Yi = a + folog 10 Xi + Si 



(4) 



where Si vanish on average, and represent the scatter of 
individual measurements about the mean relation. The 
scatter comes partly from measurement error but in most 
cases the largest contribution appears to be intrinsic (e.g. 
[IT])- It is standard practice to model the Si as being 
normally distributed with mean and some variance a 2 
i.e. Si = a5i where 5i ~ N(0, 1). We refer to this as the 
Gaussian scatter model. If the high frequency light has 
been subject to strong mixing with an ALP we expect: 



Si = aSi - log 10 C l + fi, 



(5) 



and Cj has p.d.f. /c(c) as given above. /1 is the ex- 
pectation of log 10 Ci, so the Si still have mean 0; /x can 
always be absorbed into a redefinition of the fitting pa- 
rameter a. We call this the ALP strong mixing (ALPsm) 
scatter model. The distribution of the log 10 d is both a 
distinct feature of strong mixing and very different from 
a normal distribution. Provided the variance of the in- 
trinsic Gaussian scatter (a 2 in both models) is not too 
large, it is possible, with enough measurements, to use 
the distribution of the scatter to constrain, detect or rule 
out such strong mixing. We do this by means of a like- 
lihood ratio test, comparing the null Gaussian hypothe- 
sis with the ALPsm hypothesis. Both models have the 
form Si — aSi + log 10 ((l — /) + /C;); / parametrizes the 
fraction of light that is strongly mixed. < / < 1 corre- 
sponds to partial strong mixing. However, along a given 
path, either strong mixing occurs or it does not; the X 
or 7 ray light from an object cannot be partially strongly 
mixed. / is not therefore a free parameter to be fitted. 
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The likelihood, Lf(a,b,cr,po) of the model with general 
/is: 

1 f 1 z , ? 

Lf(a,b,cr,p ) = TT / e~^/ c (c;p )dc. (6) 

( V27rcr Jo 

where Zi = log 10 F 4 -a-Mog 10 X,-h(c; f) and fo(c; /) = 
log 10 ((l — /) + fc). We fix / and po and find the values 
of a = a, 6 = 6 and cr = ct which maximize Ly. We define 
Lf(po) = Lf(a,b,a,po), and calculate 

r/(po) = 2 log ^L/(p )/-£o) • 

Keeping po fixed, both the Gaussian and ALPsm models 
have the same number of free parameters. This means 
that ri(po) is equivalent to the Bayesian Information Cri- 
teria commonly used for model selection. Conventionally, 
r i(Po) < — 6 (ro > 6) would be 'strong evidence' against 
(for) the ALP strong mixing model over the Gaussian 
one. || > 10 corresponds to 'very strong evidence'. 
If ALPs arc preferred, a useful check is to ensure that 
rf(po) is maximized for / ~ 1. If this is not the case, we 
would conclude that whilst the data is not compatible 
with simple Gaussian scatter, it is also not particularly 
indicative of the strong mixing with ALPs. 

Luminosity relations of the required form exist for 
Gamma Ray Bursts (GRBs) [TT], Blazars [H [JJ and 
AGN [JJ] and are suitable for our analysis. Addition- 
ally an 0(1) fraction of such objects are expected to be 
located within galaxy clusters. 

The 7-ray luminosity, L 7 , and radiated energy, E~, 
of Gamma Ray Bursts (GRBs) have been found to be 
correlated with a number of spectral features, giving 5 
seemingly independent relations (see [HI). Additionally 
the L 1 of Blazars, a class of AGN, is correlated with 
both their radio wave, L r |12j and near infra-red (IR) 
13J luminosities, L^. We analyzed observations of 69 
GRBs Ref. [TT], with redshifts z = 0.17-6.6, 95 EGRET 
observations of Blazars with z w 0.02 — 2.5 Ref. [T2] for 
the radio relation, and 16 Blazars (z w 0.3— 1) for the IR 
relation [13] . For all the relations however, either data 
points were too few or the intrinsic Gaussian scatter was 
too large to constrain ALP mixing. For 4 out of 5 GRB 
relations and both Blazar relations, r\ > but in all cases 
|ri| < 0.75; the sum of these ri's is only r^ RB ' Blazar as 1.6 
for po = (with similar values for other po) a statistically 
insignificant preference for the ALPsm model. 

There is also a strong correlation between the 2keV 
monochromatic X-ray luminosity, Lx, of AGN and their 
monochromatic optical luminosity, L a (at 2500A; lu m 
4.95 eV) [14]. This relation is of the form: \ogLx ~ 
a + 6 log L . We use observations of 77 optically selected 
AGNs with z = 0.061 - 2.54 from the COMBO-17 and 
ROSAT surveys as tabulated in Ref. [14 to analyse the 
scatter in this relationship. For < po < 0.4 we find: 

n(p ) AGN «l4, (7) 



and for all poj r i(Po) AGN > H- Additionally rf(p ) is 
strongly peaked at / = 1: ri(Q) w 14.1, r .95(0) w 8.6 
and r . 9 (0) « 4.9. 

There is clearly a structure in the scatter fitted better 
by ALP mixing than by the Gaussian scatter model. It 
is not clear, however, whether this is due to the success 
of the former or the failure of the latter which was only 
adopted for convenience, and because for other relations 
(e.g those of GRB or Blazars), it provides a good fit to 
the scatter. It may be that if AGN physics were better 
understood, a null hypothesis for the scatter distribution 
would be predicted that is a better fit than the ALP 
mixing model. Whilst we cannot rule out this scenario 
we can, independent of any null hypothesis, qualitatively 
check whether the structure of the scatter is really well 
matched by the ALP mixing model. 

We perform such a check by making 10 5 bootstrap re- 
samplings with replacement, D* = (x*,y*), of the origi- 
nal data, D — (xi = log 10 X%, = logi Yi). Both d and 
the D* have 77 data points. For each D* we match the 
relation s* — y* — (a* + b*x*) — by minimizing the 
0rms = h where 

fc ™(W) = (^E s * m ) > ( 8 ) 

where N p is the number of data points. k 2 is the RMS av- 
erage of the Si and is their skewness. 2d histograms 
of ki vs. kj reveal non-trivial correlations between the 
k m and, unlike the likelihood analysis, are relatively in- 
sensitive to any outlying data points. For comparison, 
we simulate data sets for both the best-fit Gaussian and 
ALPsm models in which a = 0.34 and a = 0.23 respec- 
tively. We ensure that for each original simulated data 
set ki « 0.34 and plot {fc 2 , fe} for the 10 5 resamplings. 

The detailed form of the plots varies from simulation to 
simulation, however a number of qualitative features can 
be identified as 'fingerprints' of the Gaussian or ALPsm 
models. In the former, there are two density peaks 
around {k2,k^} sa {0.34, ±0.25}, whereas in the latter 
similar peaks often occur around sa {0.23 — 0.3, ±0.15}. 
The main fingerprint of the best-fit ALPsm model, which 
occurs in most, if not all, of the simulations, is that 
most of the data points fall in a long 'tail' {k2,k^} ~ 
{0.3, -0.3} to {(0.4 - 0.5), -(0.5 - 0.7)}. These features 
persist for different values and distributions of po < 0.5, 
and in more realistic simulations where only a fraction 
(> 50%) the objects have been subject to strong ALP 
mixing. These features can be clearly seen in FIGs [TJi 
& b which are respectively typical k 2 — histograms 
for data sets simulated with the best-fit Gaussian and 
ALPsm models. Darker regions indicate higher density. 

FIG. [lp is the k 2 — k$ plot for the actual AGN data. 
There is a marked qualitative similarity between FIG.[T]; 
and the sample ALPsm model plot, FIG. [TJd. A dense 
tail-like feature is clearly present and its extent, direc- 
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Moment Plot for Simulation with Gaussian Scatter 








CO 

|-0 

o 

S -0 

"O 

E-0 

-0 
-0 

-0 




(a) 



0.4 
0.3 
0.2 
0.1 

CO 

1 o 

| -0.1 
o 

2 -0.2 
■g 

j£ -0.3 
-0.4 
-0.5 
-0.6 



0.2 0.3 0.4 0.5 
Second Moment: k2 

Moment Plot for Simulation with ALP mixing 



(b) 
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Moment Plot for Scatter in Bootstrapped AGN data 
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FIG. 1: (Colour online) Histograms of second (A^) and third 
(£13) moments (see Eq. for 10 s resamplings of 77 data points 
(a) simulated with the best-fit Gaussian scatter (a — 0.34) (b) 
simulated with the best-fit ALP mixing model (a = 0.23) and 
(c) from the observed scatter in the AGN Lx — L a relation. 
Darker regions indicate higher density. See text for discussion. 



tion and structure of peaks is typical of that seen in the 
ALP simulations. Although not shown here, there is also 
a strong similarity between the AGN k% — k$ plot and 
those found in the best-fit ALPsm model simulations. No 
evidence for any correlation between redshift and scatter 
was found, ruling out an explanation for it based on evo- 



lution of the Lx — L Q relation and / or an inaccurate 
choice of cosmological model. We note that this strong 
mixing could be independently verified by a measurement 
the polarization of 2 keV light from the AGN [5] . 

In this Letter, we have shown that the scatter in empir- 
ical X/7 ray luminosity relations can be used to constrain 
mixing between ALPs and photons. When applied to the 
AGN Lx — L relation, this shows strong evidence for 
ALPs relative to the null hypothesis of Gaussian scatter. 
Additionally, the visualizations of the AGN data reveal 
a scatter distribution with a strong qualitative similar- 
ity to that predicted by the best-fit ALP-photon strong 
mixing model. This similarity is independent of the null 
hypothesis. Strong mixing of ALPs with keV photons 
will take place in galaxy clusters if M < few x 10 11 GeV 
and 77i0 <C 10~ 12 eV, or in magnetic fields close to the 
AGN if M ~ 10 10 GeV and m < lCT 7 eV [15]. Whilst 
we cannot rule out explanations of the scatter in terms of 
known physics, it is, at the very least, a remarkable coin- 
cidence that both this and other recent analyses [5J [5] are 
fitted better by models in which very light ALPs (with 
similar couplings and masses) exist. 
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